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Abstract: We investigate the matter density fluctuations Spm/pm for two dark energy 
(DE) models in the literature in which the cosmological term A is a running parameter. 
In the first model, the running ACDM model, matter and DE exchange energy, whereas 
in the second model, the AXCDM model, the total DE and matter components are 
conserved separately. The AXCDM model was proposed as an interesting solution 
to the cosmic coincidence problem. It includes an extra dynamical component, the 
"cosmon" X, which interacts with the running A, but not with matter. In our analysis 
we make use of the current value of the linear bias parameter, 6 2 (0) = Pgg/Pmm, 
where Pmm oc (Spm/pm) 2 is the present matter power spectrum and Pgg is the galaxy 

d ' 

fluctuation power spectrum. The former can be computed within a given model, and 
the latter is found from the observed LSS data (at small z) obtained by the 2dF galaxy 
redshift survey. It is found that b\(z ~ 0) = 1 within a 10% accuracy for the standard 
ACDM model. Adopting this limit for any DE model and using a method based on 
the effective equation of state for the DE, we can set a limit on the growth of matter 
density perturbations for the running ACDM model, the solution of which is known. 
This provides a good test of the procedure, which we then apply to the AXCDM model 
in order to determine the physical region of parameter space, compatible with the LSS 
data. In this region, the AXCDM model is consistent with known observations and 
provides at the same time a viable solution to the cosmic coincidence problem. 



1. Introduction 



For the last two decades or so, the 90-years-old history of the cosmological constant 
problem Jl] has turned into the history of the dark energy (DE) problem ||. A massive 
attempt has been underway in recent times to supersede the cosmological term A in 
Einstein's equations by a variety of different entities which come under the mysterious 
name of DE, i.e., a purported new substance or, for that matter, an effective cause that 
is responsible for the observed accelerated expansion of the Universe B. In particular, 
there is the popular idea that the DE is associated with a cosmological dynamical 
scalar field (quintessence and the like) |4], which fully supplants the preeminent role 
played by A for a long time in different chapters of modern cosmology [|l], [| . Within 
this broader context, the cosmological constant (CC) has been relegated to a back 
seat in the cosmological scenario or, at least, has been degraded to represent just one 
among a host of possibilities, proposed to explain the nature of the DE. In spite of this 
situation, the A term and the standard model (ACDM model) of modern cosmology 
have been thriving rather well and have survived, essentially unscathed, the entire set 
of observational tests to which they have been exposed up to the present time (see [[J 
for a summary of the experimental situation). 

This state of affairs somehow suggests that, rather than trying to completely get 
rid of the CC term and replace it by some Ersatz entity, perhaps it would be a better 
idea to keep it and try to explain some of the unsatisfactory features of the cosmological 
standard model in terms of possible, unsuspected dynamical features of A and/or by 
introducing other dynamical complements to it. For example, while it is very hard 
to accept a small and strictly constant value of A throughout the entire history of 
the Universe, a slowly evolving DE looks more promising. Actually, this potentially 
dynamical character of the DE is the main motivation for introducing quintessence-like 
ideas . However, the contribution from the vacuum energy, most likely represented 
by the A term, is still there and remains a good candidate to be considered. Therefore, 
instead of exchanging it for a dynamical new object, it seems more economical to just 
admit that A may hide some small evolution ("running") with time or energy. Even 
more ambitiously, one may entertain the possibility that the CC could be a running 
parameter in quantum field theory (QFT) in curved space-time, i.e., an effective charge 
in the sense of the renormalization group (RG). We may call this scenario the "running 
ACDM model". An interesting proposal along these lines was developed in |7|, [8], [J 
(for a review, see e.g., ||10||), and was put to the test in [|ll|, 12[ [13| (see also JE4], |15, 16|). 



A particularly acute cosmological problem, that could be alleviated by introducing 



- 1 - 



a dynamical DE, is known as the "coincidence problem" [17], i.e., why do the matter 
and DE densities happen to be of the same order of magnitude, precisely in the current 
Universe. For example, why is the latter not one thousand or one million times smaller 
at present? The popular idea of a dynamical scalar field replacing the cosmo logical 
constant A was largely motivated by the possibility of having a framework where one 
could try to solve this conundrum. Another option for tackling this problem is to 
start with the running ACDM model, which does not deviate much from the standard 
cosmological model, and add to it a dynamical entity, X = X(t), which was called the 
"cosmon" [I8| 1 . In contrast to the standard quintessence point of view (where ACDM 
is replaced by XCDM) , here the A term is not substituted by a scalar field. Instead, it 
is assumed that there may be multiple sources of DE (perhaps some fundamental and 
others effective), including the vacuum energy, which is tied to A. The other may be 
collectively represented by the effective entity X. In this model, which was called the 
AXCDM model in ||18|| , matter and the total DE are conserved separately. However, 
the DE density here is not just pa = A/(8ttG), but is the sum of p\ and the cosmon 
density, px, i.e., Pd = Pa + Px- This total DE is locally and covariantly conserved with 
the expansion of the Universe. The advantage of upgrading the ACDM model into the 
AXCDM is that it allows for the possibility of dynamical interplay between X and A 
within a scenario where the total matter and DE densities are individually conserved 2 . 
This interplay is essential in order to provide a solution for the coincidence problem 
as well as to allow for the AXCDM model to mimic the standard ACDM model at the 
present time. Indeed, for a wide range of cosmological redshifts (including the full span 
accessible to supernovae data), the effective DE pressure and density {pd,Pd) in the 
AXCDM model may simulate a constant A behavior, pp ~ — pp, to an arbitrary high 
degree of approximation |18[ . 

It should be emphasized that, in contrast to the quintessence point of view, the 
entity X in the AXCDM model need not be a scalar field. In fact, no physical sub- 
stratum (e.g. a physical fluid) is assumed behind it. The essential condition defining 
X is the local covariant conservation of the total DE, p^ = Pa + Px, independent of 
matter. In particular, X could be the effective behavior of a more complete theory, 
comprising Einstein's gravity as a particular case. For example, it could represent the 



1 See |T^, for alternative cosmological models involving the cosmon and both variable cosmo- 
logical term A and Newton's coupling G. 

2 The name cosmon was first introduced in [plf . Here we use it in a generalized sense for any addi- 
tional component (s) of the DE, other than A, provided that the total DE density remains covariantly 
conserved. 
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net behavior of higher order terms in the effective action. 



An essential aspect to address in any proposed model for the DE is its ability to 
reproduce the existing observations. The AXCDM model has been already put to the 
test concerning nucleosynthesis bounds and the equation of state (EOS) behavior at 
small and large redhifts, that are relevant to supernovae and CMB data ||18| . However, 
one of the most important tests yet to be made is the verification of the existence of 
a region of parameter space that is also compatible with the data on the large scale 
structure (LSS) formation. This is, in fact, the main aim of the present paper. The LSS 
is contained in the galaxy fluctuation power spectrum Pgg which has been measured 



by the 2 degree Field Galaxy Redshift Survey (2dFGRS), for example, f22j. It is re- 



markable that the joint 2dFGRS and CMB analysis presented in |23[ shows that there 
is a good agreement of the LSS data with the measurements of the CMB anisotropies 



for the ACDM model as well as with numerical simulations of galaxy formation p4 [. 

On the theoretical side, this data must be reproduced by the predicted mat- 
ter power spectrum, Pmm, for any successful model of structure formation. There- 
fore, for very large scales, it is to be expected that the linear bias parameter, b 2 = 
Pgg I Pmm @, should behave as a definite, scale- independent, quantity at small cos- 
mological redshifts, i.e., when the distribution of galaxies had enough time to be cor- 
related with the mass distribution in the Universe. This is, in fact, what is obtained in 
the successful standard ACDM model, where the observed present value of b 2 turns out 



to be 1 within a 10% accuracy [22, 23]. It seems reasonable to extend this criteria to any 
other cosmological model aiming at a good description of the presently observed LSS. 
For example, in the power spectrum Pmm for the running ACDM model 0, B], which 
has been fully studied in |f?5 |, the matter fluctuations have been solved in a framework 
where they are coupled with the perturbations in the DE, in this case represented by 
the running A, described by the parameter v. Its comparison with the galaxy fluctua- 
tion power spectrum Pgg |22] puts stringent limits on the fundamental parameter v of 
this model. A non-vanishing value of v produces a time evolution of A. The explicit 
solution of the model shows that values \v\ > 10~ 4 are most likely excluded |25| . 

Here, we wish to address the LSS test for the AXCDM model. The test includes 
three parameters: v (the time evolution of A), wx (the effective barotropic index of 
the cosmon), and (the current CC density in units of the critical density). In this 
model as well, a non- vanishing value of v entails an evolving A, but in contradistinction 



to the running ACDM model or the interactive quintessence models [gg], there is no 
crosstalk between matter and DE, which is an attractive feature. The exact analysis 
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of matter density perturbations is more complicated in the AXCDM model, and we 
tackle it using the following two-step "effective method". First of all, we note that 
it is possible to ascribe an effective EOS to a variable A model. This has been ex- 
plored in detail in ^3[ . Using these results, we apply the effective EOS approach 
to the matter perturbation equations following []27|, This enables us to obtain 

an approximate treatment of the growth of matter perturbations, in which the DE 
perturbations are neglected and all the DE effects are encoded in the effective EOS, 
Pd = w e pn and in the ratio of DE to matter densities, r = Pd/pm- Secondly, to ob- 
tain useful bounds on the parameters of the model, we use the linear bias parameter, 
b 2 (z) = Pqg/Pmm- Specifically, we impose the condition ("F-test" [|30], that its 
value cannot deviate from the ACDM value by more than 10% at z — 0. From the 
above, the F-test should be essentially equivalent (although not identical) to requiring 



that 6 2 (0) = IzhO.l [22, 23|. Some concrete applications of this test can be found in |30 
In the present paper we look for the viable physical region of parameter space for the 
AXCDM model, using the aforementioned effective method. However, to check its effi- 
ciency when applied to non-trivial models with variable A, we first address the running 
ACDM model. It is important to emphasize that, in the running ACDM model, the DE 
(represented by the variable CC) and matter fluids exchange energy and, therefore, are 
interacting components of the cosmic medium. Fortunately, since we know the results 
of a full-fledged treatment of density perturbations (of both matter and DE) in this 



model |25| , we can compare them with those obtained from the effective approach, in 
which DE and matter are treated as conserved, non-interacting components and the 
perturbations of the DE are neglected. This effective approach is meaningful, provided 
we arrange that the expansion history is the same as that of the original running ACDM 



model. We may call this alternative representation the "DE picture" |T3| because it 
is close to the standard quintessence representation of the DE[j5]. The comparison 
provides an excellent test for the effective method used in the DE picture. Finally, 
we apply the effective EOS procedure, in combination with the F-test, to the more 
complicated situation of the AXCDM model and obtain the corresponding region of 
parameter space that is compatible with the LSS data. The existence of this region 
strengthens, once more, the viability and likelihood of the AXCDM model as a robust 
solution of the coincidence problem [|18|, [L^, [2(J . 



The structure of the paper is as follows. In the next section, we review the effective 
EOS approach to the computation of matter density perturbations and discuss the F- 
test. In section [| we apply the effective EOS approach and the F-test to the running 
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ACDM model. In sections |] and [5] we describe the AXCDM model and present a 
detailed numerical analysis of the matter density perturbations in order to determine 
the physical region of parameter space. The final section is devoted to discussion and 
conclusions. 



2. Dark energy, density fluctuations and the F-test 

In this section, we discuss the effective approach to the computation of the linear matter 
density fluctuations for dark energy (DE) models in which the matter components are 



canonically conserved p7| , |2q , |29f , as well as the definition of the F-test |30| . The main 
aim is to lay the groundwork for section |3], where we apply them to a cosmological 
model with a variable cosmological term A, for which an effective equation of state 



(EOS) can be defined [O, [TBI. For a general DE model within the flat Universe, the 



Friedmann equation can be written in terms of the normalized matter and DE densities 

as 

^l = Q M (a) + Q D (a), (2.1) 

where H is the present value for the Hubble parameter and a(z) = 1/(1 + z) is the 
scale factor in terms of the cosmological redshift z. The matter and DE densities are 
normalized in terms of the present critical density p° c = 3Hq/8ttG: 

M.)^, n D (. )s ^. (2.2) 

Pc Pc 

It will prove useful to also define the set of "instantaneous" cosmological parameters 
at cosmic time t, when the scale factor was a = a(t), 

ft, („) = «»(«) = (2.3) 

Pc{a) pc{a) 

where p c (a) = 3H 2 (a)/8ir G is the critical density at the same instant of cosmic time 
t. These parameters should not be confused with those in (|2.2|). We will use both sets, 
depending on the situation, and for this reason, the parameters in fl2.3| ) are written 
with a tilde. Notice that only these parameters satisfy the normalized cosmic sum rule 
at any time: 

n M (a)+n D (a) = l. (2.4) 

While the parameters of the original set ( |2.2|) also add up to one at t — to (present 
time), they satisfy the relation ( |2.1| ), in general. 
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If matter and DE are individually conserved ("self-conserved"), we can split the 
overall conservation law into two equations, 

Pm(«) + -pAi{a) = 0, (2.5) 
d 

p' D {a) + ^(l +We )p D ( a ) = 0, (2.6) 

where the primes indicate differentiation with respect to the scale factor, /' = d f /da. 
In the equations above we have parameterized the relation between the pressure and 
density of the dark energy and matter as follows: 

p D = w e p D , p M = , (2.7) 

where w e is the effective equation of state (EOS) "parameter" of the DE. It is not, in 
general, a constant parameter, but a function of the scale factor a, namely w e = w e (a). 
Therefore, we do not expect the DE to behave as a simple barotropic fluid. The solution 



of (|2.6j ) can be expressed in terms of the normalized DE density defined in fl2.2|) : 

n D (a) = Q° D exp | -3 J" ^ [1 + w e (a')]\ , (2.8) 

where Q° D = Q^(a = 1) is the normalized DE density at present. Therefore, the EOS 
coefficient is obtained from 

w e (a) = -1 - - ^ > . (2.9) 

6\lD{a) da 

Notice that we consider only non-relativistic, pressureless, matter because our pertur- 
bation calculation refers only to the epoch of structure formation. Since the matter 
conservation law is decoupled, the normalized matter density following from ( |2.5| ) reads 



Q M {a) = Q° M a- 3 , (2.10) 

where Q° M is the present normalized total matter density. 

Our analysis of matter density perturbations applies after the radiation dominated 
era. Thus, we take z (alternatively a) from about the recombination epoch z iec ~ 1100 
(i.e a ~ 10~ 3 ) to z = (a = 1, today). Moreover, it applies only to sufficiently large 
scales, where the perturbations follow the linear regime, as we discuss below. On scales 
within the horizon, fluctuations in the dark energy density disperse relativistically and 
the DE component becomes smooth, i.e., the density perturbations 5po can be con- 
sidered negligible. In these conditions, the evolution of the matter density fluctuations 
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5pM with the cosmic time t can be computed from the well-known equation of motion 
in the Newtonian approximation |32|] , 



5 M + 2H 5 M -AttGpm hi = 



(2-11) 



where, again, pressure effects (in this case pressure perturbations) are neglected in 
the matter dominated era. Here G is the Newton constant, 5m = $Pm/pm is the 
fractional matter density perturbation (density contrast) and 5m = d5 M /dt. Within 
this framework, 5m is the linear growth of density fluctuations. 

The previous equation is also valid in the presence of the DE, provided that 5po is 
indeed negligible, as assumed. In this case, H is, of course, affected by the correspond- 
ing smooth background density contribution from pp. By the same token the above 
equation is also valid in the presence of the spatial curvature term, K/a 2 , since it can 
also be treated as a smooth source (similarly to the DE). Although we will restrict 
ourselves to the flat space case, the following reformulation of ( [2.1 1| ) is valid for any 
spatial curvature. Let us write the general Friedmann's equation as 



H\a) 



8ttG 



Pm(o<) +Pd(o) 



3K 
8iiGa 2 



SttG , , 87r G Pm{o) 

p c {a) = 

3 ' 3 Q M (a) 



(2.12) 



where VLm{o) was defined in 

as 

5 M + 2H 5 M 



It is then easy to see that Eq. Q2.11| ) can be written 



H 2 Q M (a)5 M = 



(2.13) 



We wish to solve this equation for some non-trivial scenarios. First, we introduce some 
cosmetic changes in (|2.13| ), which prove very useful for practical purposes [29], as they 
allow us to apply this method to variable A models fi2| , [13|j . 

We start from ( |2.13| ), for 5m, which, interestingly enough, can be conveniently 
recast such that the effective equation of state (EOS) of the DE in the given model 
appears explicitly. We first trade the derivative with respect to the cosmic time for 
the derivative with respect to the cosmic factor: 5m = a H5' M , where 5' M = (15m/ da. 
Similarly, 5m = (a H 2 + a H) 5' M + (a H) 2 5" M . We observe that we can eliminate H (i.e. 
the time variation of the Hubble function) from the last equation by noting that, in 
the matter dominated epoch, it can conveniently be written as (henceforth, we confine 
ourselves to the flat case only) 



H = —4 7r G [ P M + (1 + w e ) p D ] = -- H 2 



1 + 



w P r 



1 + r 



(2.14) 
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where w e = Pd/pd is the aforementioned effective EOS parameter of the DE and 
r = r(o) is the ratio between the DE and matter densities at any given time |]18||: 



fiiif(a) n tf (a) Pa/ (a) 

The relation ([2.15 ) will play a relevant role in the study of the cosmic coincidence 
problem within the AXCDM model (see section |J). Substituting the previous equations 
into Eq. ( |2.13|) , we finally obtain 



KM) + - 



w e (a) r(a) 
1 + r(a) 



5' M (a) 3 1 5 M (a) 



2 1 + r(a) a 2 



0. (2.16) 



This will be our master equation to evaluate the "effective" growth of linear density 
fluctuations. As we said above, we assume that the DE component becomes smooth 
within the horizon. The DE effects enter our calculations only through the effective 
EOS function ( |2.9| ) and the ratio ( |2.15| ), and, thus, the linear growth of matter fluc- 
tuations is computed in an effective way. We also assume that the DE density was 
negligible at the recombination era and remained so until z < 10, when all relevant 
modes for LSS formation had already entered the horizon. In particular, the pertur- 
bation amplitude at recombination, when the CMB was formed, is independent of the 
particular DE model since po <C Pm at that epoch. 

To better assess the meaning of Eq. ( 2.16| ), let us first consider two simple examples. 
In the absence of DE (tip = and, thus, r = 0) the growing mode solution of the 
above equation in the matter dominated epoch (a ~ is very simple and well- 
known, 5m ~ cl. However, in the presence of DE, the growing mode solution is more 
complicated. Assuming a time interval not very large such that tip and remain 
approximately constant, we can take w e = — 1 and, hence, ( |2.16|) reads 



«,(«) + 5(1 + - |(1 - H B ) S -ZgL = 0. (2.1 7, 



Looking for a power-law solution, 5 ~ a p , in the limit Qd <^ 1? we find 



6 a 
5 



hi ~ a 1 - 6 ^/ 5 ~ a (1 - - Q D In a) . (2.18) 



This equation conveys, very clearly, the physical idea of growth suppression when a 
(positive) DE density is present within the horizon. Although tip is not constant in 
general, this qualitative feature should persist in more realistic situations. Furthermore, 
if the DE has some smooth dynamical behavior, we cannot exclude that Qd could have 
been negative for some period in the past. If so, the previous equation also shows 
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that structure formation was reinforced during that period. In this paper, we wish to 
study the exact numerical solutions of ( |2.16| ) in some non-trivial scenarios, where the 
cosmological term is not only arbitrary and non-vanishing, but is evolving smoothly 
with time. Specifically, we wish to solve ( |2.16|) for both the running ACDM and 
AXCDM models, mentioned in the introduction. The clue to doing this is to mimic 
these variable A models through a non-trivial effective EOS, w e = w e (a). 

From the above discussion, it follows that the linear growth should essentially 
behave as 5 m = a near z = z rec . We are going to use this property as an exact boundary 
condition for solving (|2.16| ), together with 5' M = 1 at recombination. Equivalently, if 
we introduce the standard linear growth suppression factor, G(a) = 5m (a) /a, Eq. ( |2.16| ) 
is readily seen to transform into 



G"(a) + 



7 3w e (a)r(a) 
2 ~ 2 l + r(a) 



g\a) + 3 [1 - w e {a)]r{a) g(o) = Q ^ ^ 



a 2 1 + r(a) a 2 

which we solve with the boundary conditions Q(z rec ) = 1 and Q\z rec ) = 0. 

To compare our predictions with observations, it is useful to invoke the linear 
bias factor. Its present value is defined as 6 2 (0) = PqgI Pmm, where Pmm is the 
computed matter power spectrum at the present epoch, starting from the Pmm at 
the recombination epoch (obtained from the observed CMB) and Pqg is the galaxy 
distribution power spectrum (obtained from the observed galaxy-galaxy correlation 
function). Therefore, b 2 (0) measures the difference in clustering between galaxies and 
mass fluctuations, i.e., it parameterizes the degree by which light traces mass. It can be 
related to the rms mass fluctuations on random spheres of radius R h~ x Mpc, typically 



with R = 8, and hence connected with erg |32||. In general, the bias factor can also 



be defined in the non-linear regime. Here, however, we are concerned only with the 



linear bias factor [23], which measures the difference in clustering between galaxies and 
mass fluctuations at very large scales, namely at scales for which the wave-numbers of 
the Fourier modes are in the range 0.02 < k < 0.15 h Mpc -1 . The observational data 



concerning the linear regime do, in fact, lie in this range |22|| . Expressing the scales in 
terms of the Hubble radius Hq 1 ~ 3000 hr 1 Mpc {h ~ 0.7), this implies values of k up 
to 450 -f^o- Thus, the minimum length scale explored by the available LSS data is of 
order of 10 Mpc ~ 8 h' 1 Mpc spheres. 

For the computation of the linear bias factor, we need the matter power spectrum, 
whose general structure is 

P MM (k, a) oc f 6 -^-) = 5 2 M (k, a) = a 2 G 2 (k, a) , (2.20) 
V Pm J 
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where <5m(&, a) and Q(k, a) are the Fourier transforms of the corresponding quantities. 
In general, the correlation function for the mass distribution need not coincide with 
the correlation function for galaxies. Rather, a "bias" between the two is expected ||32|| . 
However, at the LSS level, it is also expected that the value of the linear bias should 
be some scale-independent number in the late epochs of structure formation (z ~ 
0), namely when galaxies have had time enough to be correlated with the mass, or 
equivalently, when the gravitational pull has drifted them to overdense regions. As 
a matter of fact, the observed galaxy power spectrum Pgg-, emerging from the final 
2dFGRS catalog, did indicate these features very clearly. Most remarkably, the data 
pointed to the value b\(z ~ 0) = 1.0 to within a 10% accuracy for the ACDM model 
[pm , i.e., when Pmm in ( |2.20| ) is computed from the growth factor for the standard 
cosmological model, characterized by strictly constant A. This is in agreement with 
the previous result, b(L s ,z ~ 0) = 1.10 ± 0.08, of the 2dFGRS collaboration for the 
APM-selected massive galaxies (Lg = 1.9L*), averaged over all types [^3|], indicating 
that there is one galaxy per dark matter halo of mass ~ 10 13 h^ 1 M & at the present 
epoch. 

As indicated before, in the evaluation of the bias factor, P GG is fixed by the LSS 
data, whereas Pmm is a theoretical quantity - and, hence, model dependent. However, 
the very good prediction (at the 10% level) of the bias factor by the ACDM near 
b 2 (0) = 1, suggests the following strategy to put limits on new models of structure 
formation. Rather than computing Pmm in detail for the given model, it can just be 
compared to the ACDM model. This is, of course, simpler than the full computation 
of Pmm, in which details of the transfer function and other normalization factors must 
be included in the prefactor of (|2.20|) . These prefactors are common and cancel in 
the ratio. Therefore, in the present analysis, we adhere to the following "F-test" (first 
proposed in[[!(J), which uses the method of comparison to evaluate the viability of a 
given DE model. To gauge the deviation of the power spectrum of the model, Pmm (a), 
with respect to that of the ACDM model, P MM (a), we define the parameter: 



\F\ = 



PmmM — PMWfi. a 



Pmm{°) 



a=l 



a=l 



where Q{a) is the solution of fl2.19|) for the given DE model with some effective EOS 
w e = w e (a) (in some cases w e can be constant, but in general, it is a function of the 
scale factor), and Q\(a) is the corresponding solution for w e = —1 (ACDM). From its 
definition, the factor F is a number, computed for a = 1 (i.e., the current Universe at 
z = 0). On the other hand, since we compare all models to the same observed galaxy 
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power spectrum Pgg-, it follows that the above F-factor can be directly related to the 
relative difference between the linear bias factor of the given model with that of the 
ACDM model: 



1 - 



2 = 



b\z) 



(2.22) 



2 = 



Pmm{ z ) 

As emphasized above, observations ^3] show a scale invariant linear evolution of 
b\(z) for the ACDM model towards b\(0) = 1.0 ± 0.1 at present. Therefore, since all 
the DE models should, presumably, approach the predictions of the successful ACDM 
model near the present time, we require that any given DE model with a growth factor 
g(a) should pass the following "F-test" pi: 



1 - 



Q 2 M 



G 2 (a) 



< 0.1. 



(2.23) 



a=l 



In performing the test, both G(a) and Q\{a) must be evolved from the recombination 
epoch, where the initial conditions are fixed (see above), to the present time. It is 
understood that the maximum limit, -F max = 0.1, is saturated when Q\{1)/Q 2 {1) = 1.1 
and that the minimum limit, F m[n = —0.1, when Q\{1)/Q 2 {1) = 0.9. In the next 
sections, we apply these limits to the running ACDM and AXCDM models in order to 
bound their respective parameter spaces. 



3. Effective method approach to the running ACDM model 

The F-test, defined in the previous section, is not exactly equivalent to requiring that 
6 2 (0) = 1 ± 0.1 for a given model, but it is not very different from it and has the 
advantage of being a relatively economical procedure. However, we need to check its 
efficiency in some non-trivial situation before applying it to a complex DE model, 
such as the AXCDM model. To this end, we first apply the effective approach to the 
computation of matter fluctuations, together with the F-test, to the running ACDM 
model, for which we know the results of a complete analysis P5fl . 

In the running ACDM model, the CC, or equivalently, its associated energy density 
Pa = A/87T G, is an evolving parameter. It "runs" because of the quantum loop effects of 
the high energy fields and, therefore, it satisfies a renormalization group (RG) equation. 
We refer the reader to the original literature for details [[j], |8|, P] . Here, we just 
highlight the basic concepts and equations, needed for our analysis. In this framework, 
the physical RG running energy scale is identified with H (the Hubble parameter) and 
the solution of the aforementioned RG equation reads 

p A (H, v) = pi + — Ml {H 2 (a; v) - H 2 ) , (3.1) 
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i.e., the DE density (in this case a running CC term) is an affine quadratic law in H, 
where p\ and H are the current values of these parameters. In this model we have a 
single (dimensionless) new parameter u, given by 

M 2 



a 

v= 3.2 

12 Ti Mp ' 



which is essentially the ratio (squared) of the effective mass M of the high energy fields 
to the Planck mass (Mp), with fermions contributing a = — 1 and bosons, a = +1. If 
the effective mass M of the heavy particles, associated to some Grand Unified Theory, 
is just the Planck mass Mp, the parameter v takes the value (positive or negative, 
depending on a) 

»o = -J- ^ 2.6 x 10" 2 , (3.3) 

sometimes referred to as the "canonical" value. In practice, the preferred values of v 
are smaller than v§, as we will see, which is, in fact, the natural situation because it 
corresponds to having the heavy particles at some scale nearby, but below, the Planck 
scale. We point out that the Hubble function of the running model is also ^-dependent, 
H = H(a; v). In the flat case, 

a- 3 ^) -1" 



E\a ] v)=Hl 



1 + 0° 



1 — V 



(3.4) 



For v = 0, we recover the standard form corresponding to a strictly constant A. 

In the framework of the running ACDM model, there is energy exchange between 
the vacuum and matter sectors and we have the following mixed conservation law: 

3 

p' A (a) + p' M (a) + - p M = . (3.5) 
a 

From (|3.1| ) and (|3.5| ) we see that a non- vanishing value of v causes, not only a running 
of the CC density as a function of the scale factor (or the redshift), p\ = p\(a;i>), 
but also an exchange of energy between p A and pu- For v = 0, however, p A becomes 
constant and (|3.5|) boils down to the old matter conservation law (|2.5| ). These features 
are also apparent from ( [3TTD and (|3.4j ). 

For all the new dynamical features that a variable CC term may entail, it should 
be clear that its EOS parameter still remains w\ = —1. In this context, we may speak 
of the model as being described within the "CC picture" (TB| , that is to say, the original 
formulation, in which the A term is explicit and the matter density is non-conserved. 
However, it may be advantageous to perform a "change of picture", i.e., a description 
of the running ACDM model within the "DE picture" | I3]| . In the latter, we envision 
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the given running CC model as if it were a DE model with the same expansion history 
(that is, with the same numerical values of H) but with self-conserved matter and DE 
densities ( j2.5|) and (|2.6| ). The numerical matching of the Hubble functions is essential 
in order to guarantee that the physical results are the same in both pictures. The 
basic reason for moving into the DE picture is because we wish to find a representation 
of our cosmological system where we can compute the matter perturbation equations 
through Eq. (|2.16|) , in which the DE effects are confined only to the effective EOS 
w e = w e (a) and the ratio r(a) = p£>(a) / ' Pm{o)- Since, however, (|2.16|) stems from the 
original Newtonian form ( [2.11| ), which is derived under the hypothesis of total matter 



conservation [ 3"2| , we must use a formulation of our cosmological scenario in which this 
condition is also satisfied. This representation is provided by the DE picture. 

While this alternate formulation is perfectly possible, the fact that the matter 
density is non-conserved in the original CC picture, suggests that the mapping of 
the latter into the DE picture can only be carried out at the expense of finding a 
non-trivial effective EOS parameter (actually some complicated function w e = w e (a)) 
relating the self-conserved pr> density and pressure pa in the DE picture, pn = w e pD- 
This non-trivial EOS function for the running ACDM model was determined in fl2|, |I3||. 
For the present purposes, it is convenient to first determine the normalized DE density 
VLd = &d(z] v) explicitly and then apply (|2T9|). Following the procedure indicated above 
and adjusting the two pictures such that the current values of Qd and fl\ coincide, we 
find 3 

Q D ( a] v)=QI + 3l_ a -3 ( a 3, _ 1+u \_^3l (3 . 6) 
1 — v 1 — V 

where = 1 — Q° M . (A generalization of this method allowing for the mapping of 

any cosmological model with variable cosmological parameters from the original "CC 

picture" into the "DE picture" was developed in [O], to which we refer the reader for 

details on the entire procedure.) From this expression and (|2.9|) , we find the desired 

result, 

a , (a , v) = _ t + (1 _ ,) . _ _ ^ _ _i JL _ _ _ ^ . (3.7) 

This effective EOS behaves nicely: w e — > — 1 for a — > 1. Expanding linearly the 
previous formula for small u near our epoch, we have 

Q° A r In a . 
w e {a-v) ~ -1+3^^ — g-. (3.8) 



3 For comparison, the corresponding normalized DE density in the original CC picture is given by 



J]), where H has the explicit form (3.4), divided by p° c . See Refs. [|TlJ and [[13[ for numerical plots 
of these functions. 
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Figure 1: The effective EOS for the running ACDM model as a function of the cosmological 
redshift z. We assume flat space geometry and the prior tt° M = 0.3. The shaded area satisfies 
the "F-test" condition \F\ < 0.1, see Eq. ( ^23|) . 

Since a(t) < 1 for any look-back time, this equation shows that the self-conserved 
density pd in the DE picture for this model does not necessarily behave like quintessence 
(w e > —1) in our past; it is only so if v < 0. It, however, behaves like "phantom 
DE" (i.e. w e < -1) when u > 0. 

These features can be clearly identified in Fig. [l], where we show the evolution of 
the effective EOS for the running ACDM model, assuming a prior Q° M = 0.3 for the 
normalized matter density at present. For convenience, we display w e = w e (z) as a 
function of the cosmological redshift z = (1 — a) /a. The numerical examples shown in 
this figure correspond to four values of z/, expressed as small fractions of the canonical 
value (|3~3] ). Specifically we plot (|3.7| ) for v = ±0.1 v G and v = ±0.05 z/ . Although 
not shown in the figure, from the analytic expression (|3.7|) , it follows that for v > 0, 
the effective EOS for the running ACDM model goes to zero at very high redshift, 
whereas if v < 0, w e (z) — > —v for z — > oo. We note that if v > 0, the denominator 
of Eq.( [3.7|) vanishes at some redshift Z\ > 0. This is related to the fact that for a 
positive v, the normalized DE density (|3.6|) vanishes at some point in the past. This is 
easily seen from the sign of this function, which changes from Qd{(1', v) < in the far 
past (Q D — > —Q° M /a 3 < 0, for a <C 1) to 0, D (a; u) > at present (Qd > 0, for 

a — > 1). This transition is not possible for v < 0. The divergence of w e for positive v 
is not a real singularity of the theory since the fundamental cosmological functions in 
both pictures, pu = Pm(cl), pa = pa(o),Pd = Pd(cl), H = H(a), etc., are finite for all 
a. The discontinuity of w e = w e (a) is associated with the effective description of the 
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Figure 2: F-test for the running ACDM model under the same conditions as in Fig. [T]. 
We show the factor F, Eq. ( 2.2l| ), versus v. This parameter becomes confined in the range 
\v\ < 0.05 z^o — 10~ 3 , due to the condition \F\ < 0.1, defined in ( |2.23| ). The horizontal dashed 
line corresponds to the maximum allowed deviation, \F\ = 0.1, of the bias factor of this 
running model with respect to that of the standard ACDM model. 



original model in the DE picture and has no effect on our analysis. Indeed, the product 
function w e (a) r(a) in the differential equations Q2.16Q and fl2.19|) remains finite for all 
a. Therefore, the computation of the growth factor is free from singularities in the 
entire range of definition. 

Following the procedure explained in detail in section |2|, we have determined the 
growth factor for both the standard ACDM and running ACDM models. For the 
former, we naturally used w e = —1, and for the latter, Eq. (|3.7|) . After solving the 
perturbation equation ( p.!9| ) numerically, we computed the parameter F, defined in 
(p.21|) , and applied the F-test ([2.230 , estimating the range of allowed values for the 
parameter u, shown in Fig. |2|. The two values of v which saturate the F-test are the 
following: the upper bound, -F max = 0.1, which corresponds to z/ max ~ —1.3 x 1CT 3 and 
the lower bound, F min = —0.1, which corresponds to u ui \ n ~ +1.4 x 10 -3 . They are 
represented by the two black circles in Fig. |^. Therefore, the parameter v is confined 
to the range \v\ < 0.05 u = 10~ 3 , corresponding to the shaded areas in Figs. [I] and 
|2|. The other values, v = +0.1v and v = — 0.1u , considered in Fig. [I], correspond to 
F min ~ —0.17 and F max ~ 0.20, respectively, and are excluded because they are out of 
the allowed limits (the shaded region). 

According to the 2dFGRS collaboration, the central values of the normalized mat- 
ter and DE densities at present, obtained from the combined LSS and CMB data for 
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the flat ACDM Universe, are Q° M = 0.27 and Q° A = 0.73 respectively |2]§. These re- 
sults are in good agreement with the most recent analysis by the WMAP team, which 
reported Q° M = 0.28 and Q° A = 0.72 |J. To calculate the limits on v in this work, 
we considered flat space with the following values of the matter and DE densities at 
present: fl° M = 0.3 and Q° D = 0.7, respectively. 

We end this section by noting that the limits we obtained for the fundamental 
parameter of the running ACDM model, 

-10~ 3 < ^<+10~ 3 , (3.9) 
are essentially in agreement with previous estimates obtained by a number of different 



methods and authors. In particular, the results of |yj are in good agreement with (|3.9|) , 
although in the latter study an indirect procedure was employed, which was based on 
bounding the amplification of the matter density spectrum in the recombination era 
caused by vacuum decay into CDM. We emphasize that our results are also in good 
agreement with the full calculation of coupled perturbations of matter, metric, and 
DE for the running ACDM model, presented in [p5|] , where the DE fluctuations were 
included as perturbations in the A parameter. A tighter bound for v was obtained, 
\v\ < 10~ 4 . However, taking into account the economy of our present procedure, as 
compared to [p5|| , we can assert that the effective approach to the calculation of the 
matter density fluctuations, in combination with the F-test, provides a reliable estimate 
of the physical region of the parameter space. 

From the tiny range of values of v obtained, ( |3.9|) , it is patent that a procedure 
based solely on the search for non-standard features in the EOS profile (i.e. departures 
from the featureless behavior w e = —1) can only be efficient if sufficiently large redshift 
observations are used, which may be able to be obtained from the DES and SNAP 



programs ||35|| . Even these programs may have difficulty in differentiating the running 



AXCDM model from the standard ACDM model, where w e = —1.0. From Fig. |l], 
we see that the allowed interval for the effective EOS parameter of this model at the 
redshift z = 1.7 (the largest one reachable by SNAP) is rather small, and w e does not 
deviate appreciably from -1.0: 

-1.03 <w e {z = 1.7) < -0.97. (3.10) 

The effective EOS parameter for the running ACDM model varies very slowly with z. 
For example, at z — 2, —1.05 < w e (z — 2) < —0.96. If a small departure from the 
strictly constant w e = — 1 could be observationally substantiated, it would be a strong 
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sign of dynamical behavior of the DE. A remarkable feature of the running ACDM 
model is that its effective EOS can mimic departures from the standard cosmological 
model, both in the quintessence and phantom regimes, even though the original model 
has nothing at all to do with quintessence or with phantom-like fields. The original 
model (in the "CC picture") is nothing more or less than a model with a true, albeit 
running, cosmological term and, hence, with wa — — 1- 

The most important result of our investigations of the running ACDM model is 
that the LSS data are strongly sensitive to the running of the A term, which is, of 
course, the reason why the parameter v becomes strongly constrained. This confirms 
the results of [^5| within the present effective approach. However, as noted earlier, 
in this model, the variability of A is entangled with the non-conservation of matter. 
Therefore, in the next section we investigate a model (the AXCDM model) with a 
dynamical A where matter is conserved, to see if the limits on v can be more relaxed. 
In the light of the successful application of the effective method, we apply the same 
approach to the AXCDM model. 



4. The framework of the AXCDM model 

In the introduction, we have already explained the basic motivation for the investiga- 
tion of the AXCDM model, i.e., the possibility of its being a solution to the cosmic 
coincidence problem 4 . To achieve this important property, the DE density in this 
model must have two components, to wit: the running p\{a\v), which satisfies the 
RG evolution equation ( |3.1| ), and the dynamical component X (the "cosmon"), whose 
underlying nature can be very general and need not be specified here. These two 
components of the DE may be in interaction. As already stated, X is not necessarily 
related to any scalar field. It could be some effective entity, related to the structure 
of a more general theory, in which Einstein's gravity is embedded, and which includes 
the effect of higher order terms. It is supposed that the X component behaves as a 
barotropic fluid with a constant EOS parameter wx = Px/px- (In general, wx may 
not be a constant, but a function of redshift.) Typically, the index wx for the cosmon 
is in one of the following two expected ranges: wx ^ — 1 (quintessence-like cosmon) or 
Wx < ~ 1 (phantom-like cosmon) 5 . Adopting the simplest possible AXCDM scenario, 



4 We refer the reader to a detailed discussion of this model in |l8|, Q ■ 

5 As noted in [jl8| , quintessence-like and phantom-like cosmons do not necessarily exhibit the naively 
expected behaviors corresponding to quintessence (dpx/da < 0) and phantom energy (dpx/da > 0), 



respectively, since in general there is an interaction between X and the running A, see (4.2) 
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we assume that the total matter density does not interact with the DE and that it is 
covariantly conserved, thus satisfying ( |2.5| ). If we define the total DE density as the 
sum of the CC density and the cosmon density, 

p D (a) = p A (a) + p x (a) , (4.1) 

it follows from our assumption of matter conservation that pn is also covariantly con- 
served. The cosmon X is actually defined through this conservation condition [TE[ . 

Therefore, the quantity (|4.1|) satisfies Eq. fl2.6|) , which can be recast as 

3 

PM + p'x(a) + -(1 + w x ) Px{a) = . (4.2) 
a 

From ( ^4.2| ), it is clear that, although the total DE density is locally conserved in the 
AXCDM model, in general, the individual px and pa densities are not. There is a 
transfer of energy between them, which is governed by the above equation. 

The expansion history of the AXCDM model is determined by its Hubble function 

H\a) = Hi [n° M a' 3 + tt D (a)] , (4.3) 

where Ho(a) = p£>(a)/ p° c is the normalized total DE density ( |4.1| ). The corresponding 
expression that satisfies the above equations in the matter-dominated, flat Universe is 
given by 



1 — v Wx — e 



i - n\ Wx n° M 



1 — 1/ Wx 



a -3(l +wx -e) ^ ( 4 4 ) 



where, for convenience, we have defined 

e = u(l + w x ). (4.5) 

We will see, below, that this quantity must remain small for the model to be compatible 
with primordial nucleosynthesis. 

The normalized densities, at present, satisfy the relation ( |2.1| ). For a = 1 (i.e. 
z = 0), it takes the form 

n D + n° M = n o x + nl + n° M = i, (4.6) 

which is the current cosmic sum rule. With the help of this relation, it is easy to see 
that, for v = 0, the DE density ( f4.4j ) boils down to 

n D ( a ) = n° A + n° x a - 3il+wx) . (4.7) 

Clearly, in this particular case, where u=0, the AXCDM model mimics a system, 
consisting of a quintessence or phantom fluid (depending on the value of wx), together 
with a constant cosmological term. From ( |4.7p , we then have three possible scenarios: 
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i) If Wx ^ — 1 and Q x > (quintessence-like cosmon), the expansion of the 
Universe can be stopped, provided that f2° < 0, since the X density gradually 
diminishes with the cosmic time and the r.h.s. of (|4.7|) becomes negative. Hence, 
there exists a future time, a = a* > 1, when if (a*) = 0. As mentioned previ- 
ously, in the AXCDM model, the total Q° D ~ 0.7 (corresponding to Q° M ~ 0.3). 
Therefore, due to the contribution from Q x in the sum rule ( |4.6| ), the possibility 
that < should not be discarded a priori; 

ii) If wx ^ — 1 and Q x > (phantom-like cosmon), the DE density increases 
indefinitely with the expansion. It does not matter whether Q A is positive or 
negative, the Universe unavoidably ends up in a super-accelerated phase, which 
triggers a catastrophic disruption of all bound systems, a singularity known as 
the "Big Rip" p 



iii) If Wx < — 1 and Q° x < 0, the cosmon density, although phantom-like, acts 
with matter to decelerate the expansion of the Universe. We have the opposite 
situation to the Big Rip: the Universe becomes super-decelerated. The kind of 



cosmon able to create this scenario was previously called "phantom matter" [18 
Phantom matter, therefore, avoids the Big Rip and helps the Universe to reach 
a = a* > 1, where if(a*) = (i.e., a stopping point). In the present instance, 
this point will exist provided > 0. Obviously phantom matter is special 
in that it corresponds to negative energy density, which is, however, not new 
in the literature |34[] . In spite of its rather peculiar nature, phantom matter 
satisfies the strong energy condition (see Fig. 1 of |18|). As previously discussed, 



the cosmon may well be an effective entity and, therefore, could simulate the 
behavior of phantom matter. This is in contrast to the "standard phantom 
energy", considered in the previous case, which violates all of the classical energy 
conditions and leads to a cosmic doomsday. 

From the previous examples, with v = 0, it is clear that there are simple scenarios 
within the AXCDM model, in which the cosmological expansion can be stopped at 
some point in the future. The "stopping" point is actually a "turning point" in the 
evolution of the Universe; it bounces back at that point and is subsequently redirected 
towards the Big Crunch. However, stopping can be formulated on very general grounds 
within the parameter space of the AXCDM model and is not restricted to v = 0, as in 
the previous examples. This issue is central to the cosmic coincidence problem |Tj| [HJ 
and is correlated with the existence of a maximum of the ratio r(a), defined in ( |2.15| ), 
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which gives the amount of DE energy versus matter at any time. To further address 
this problem, let us compute explicitly the function r(a) in the AXCDM model in the 
matter dominated era. With the help of (|4.4j ) and (|2.10|) , we find 



r(a) 



- u) a 3 



1 



+ 



lw x 



+ 



0° 



w x 



^(1 



-3(w x -e) 



(4.8) 



v ) Wx — e 

In order to provide an acceptable explanation for the cosmic coincidence problem, this 
ratio should be bounded and stay relatively small throughout the entire history of the 
Universe. This can be expressed through the condition 



r(a) 
r 



< 10 



(4.9) 



where ro = Q° D /Q° AI ~ 7/3 is the present value of r (of order one). The ratio (fO 



should, therefore, reach a finite maximum in the evolution of the Universe. The con- 



ditions for this to occur can be expressed as |j 

> 0, 



0° 
"a 



V 



[i + wx) {n° A -v) <o 



(4.10) 



w x (n x + vn° M )-e(i-ni) 

We can easily check that the simple stopping scenarios i) and iii), mentioned above, 
are consistent with these requirements. We note that wx = — 1 (i-e. a pure CC-like 
cosmon) is not allowed. 

Besides the two conditions ( fO| ) and ( |4.10| ), the ratio r(a) should satisfy the nu- 
cleosynthesis constraint, namely that its value at the primordial nucleosynthesis epoch 
should be |rjv| < 0.1 [^7fl . In that early epoch, the ratio r(a) is no longer given by (fi~8l ) 
since we must use the radiation equation for the (relativistic) matter density, namely 
^i?( a ) = ^?j a ~ 4 instead of ( p.lO| ). Thus, we have: 



w R -w x + e 



+ 



1-^a 
0^(1 - u) 



Wr - W X 



(I 



N 



(4.11) 



w R -w x + e 

where w R = 1/3 is the barotropic index for radiation and ax ~ 10 -9 , the scale factor 
at the nucleosynthesis time. Since wx < —1/3, the condition 



lei < 0.1 



(4.12) 



insures that the contribution from the second term on the r.h.s. of ( [4.11| ) is negligible 
and we find that ( |4.12|) is essentially equivalent to |rjv| < 0.1. Hereafter, we shall refer 
to (|4.12|) as the nucleosynthesis bound. 



For v 



and n° x 



0(O° 



reduces to 



r(a) 



0° 

_^A_ 3 



(4.13) 
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which is, of course, the standard ACDM model prediction for the ratio r(a). We see 
that, as the time passes, a —>■ oo and, therefore, this ratio is unbounded from above, 
i.e., it can take any arbitrarily large value. Thus, the cosmic coincidence problem boils 
down to understanding why, at a = 1 (the present time), the ratio is just of order one. 
In other words, what makes our time special? This question has no reasonable answer 
within the standard ACDM model. It also has no acceptable answer within the running 
ACDM model. Moreover, the ratio r(a) in the running ACDM model cannot be worked 
out as a particular case of the AXCDM model because matter is not conserved in the 
former, whereas it is conserved in the latter. 

Let us further elaborate on this ratio by considering its evaluation within the 
running ACDM for the two cosmological pictures that we are considering, namely the 
CC and the DE pictures. In the former, we have r(a) = p\(a)/ 'pm{°)i where Pm{o) 
and p\ (a) can be obtained from equations (7)-(9) of Ref. ||, for example. The final 
result reads 

r(a)cc = |_ ^ + JL. [, _ a 3 M] = „3 M + _J_. (4 , 4) 

For Vt\ < v, the CC density eventually becomes negative and, in that case, there is 
stopping in the running ACDM. However, since f2° ~ 0.7, this possibility entails a 
value of v, which is ruled out by our result (|3.9|) . Hence, there is no viable solution 
to the cosmic coincidence problem in this model. We wish to stress that this result 
is independent of the particular cosmological picture chosen to derive it. Indeed, in 
the DE picture the ratio r(a) is, instead, r(a) = Pd{o)/ Pm{o) = &d{p)/Qm{q), where 
VLm{.q) and VLd{o) are given by (|2.10| ) and ( |3.6|) respectively. Therefore, 



Q° A v \ „ g? v - 1 + v Sl\-v 



a 



3v 



where we used VlP M = 1 — . Notice that both ratios, r(a)DE and r(a)co ; satisfy the 
correct normalization at the present time, i.e., r(a = 1) = Q^/Q^j. Again, we see from 
( j4.15|) that the condition for the DE density to become increasingly negative (leading 
to stopping) will occur only if Q° < v. Thus, we obtain the same conclusion as in the 
CC picture. 

The foregoing results indeed show that irrespective of the cosmological picture 
used to perform the analysis, the conclusion is the same, to wit: in the running ACDM 
model there is no natural solution to the cosmic coincidence problem. Although this 
model does provide some interesting dynamics for the CC term, it does not have the 
ability to ameliorate the cosmic coincidence problem. It is only when the X entity is 
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introduced in interplay with a dynamical A, that the ratio r(a) takes the form of ( O ) 
and can be kept within bounds throughout the entire history of the Universe, which is 
indeed what is needed to solve the cosmic coincidence problem. 

The three conditions, flOD , flllOD and ( |Ll2|) , define a limited 3D -region in the 
parameter space ($7°, Wx, v) of the AXCDM model at a fixed Q° RI ~ 0.3 [18|. In 



the next section, we give a detailed numerical analysis of the AXCDM model, after 
including two more conditions, both of which are related to the effective EOS of the 
model. One of the conditions is its compatibility with the LSS data. We express this 
compatibility condition, using the F-test ( |2.23|) on the linear bias parameter of the 
model, which depends on the effective EOS approach to the calculation of the growth 
factor. The other condition is the maximum allowed deviation that we can tolerate 
for the value of the effective EOS parameter, w e = w e (z), away from the CC value, 
w e = —1, at z — 0. As we shall see, with these five conditions, we will be able to 
significantly improve the determination of the physical volume in the 3D parameter 



space of the AXCDM model, as compared to [18 



5. Numerical analysis of the AXCDM model: cosmic matter 
perturbations versus cosmic coincidence 

In this section, we present a complete numerical analysis of the AXCDM model, in 
which its most salient feature is the inclusion of the effective EOS approach to the 
growth of matter density fluctuations, together with its ability to solve (or at least sig- 
nificantly alleviate) the problem of cosmic coincidence, mentioned in the introduction. 



The effective EOS for this model is [18 



, , Pa + Px -Pa + wxPx -i , , v fijr(a) 
w e {a) = — = = -l + (l+u> x )— — - (5.1, 

Pa + Px Pa + Px HdW 
where = fi D (a) is given by (|4.4|), and 



n x (a) =[n x STL a - 6(l+wx ~ t} + n u M a-' . (5.2) 

V w x -e J w x - e 

Note that for v = , Eq. (]5.2|) reduces to the expression in the second term on the r.h.s. 
of ( |4.7p , as it should. Equations Q5.2Q and ( |4.4j) may be used to compute the effective 
EOS ( |5.1|) . Alternatively, (|4.4]) may be substituted in ( p.9|) , the two procedures being 
equivalent since 

1 + w x ) fix a = -7T ■ 5.3 

3 da 
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Clearly, the dynamical features of the cosmon play a preeminent role in the behavior of 
the effective EOS for the AXCDM model. For example, if the cosmon is quintessence- 
like (wx ^ — 1 and f^-(a) > 0) near a = 1, the overall effective EOS of the model 
will be quintessence-like near our time. However, if the cosmon behaves as phantom 
matter (viz. wx ^ — 1 and Qx{ a ) < 0), the effective EOS of the model will still be 
quintessence-like. On the other hand, if the cosmon behaves as standard phantom DE 
(wx < —1 and fix (a) > 0), the AXCDM model will also behave phantom-like as a 
whole. In general, the AXCDM model behaves effectively as quintessence (phantom) 
near z = 0, if and only if 1 + wx and Qx have the same (opposite) signs. This can be 
seen from fl5T3| ) and (|5.1[). 

For the numerical analysis, we insert the EOS formula ( [5.1|) into ( |2.19|) for the 
effective growth of the density perturbations. We then solve this equation and perform 



the F-test ( 2.23 ). That is, we impose a condition on structure formation, whereby we 
discard all points of the parameter space, for which the growth factor of our model 
at z = deviates by more than 10% from that of the ACDM. In the analysis, we 
have to also include the three conditions discussed at the end of the last section, 
namely, the primordial nucleosynthesis constraint, the stopping condition, and the 
bounding condition on the ratio r = r(a). We begin the analysis by checking the 
ability of the F-test, alone, ( |2.23| ), to define a limited region of the (wx, v) plane. As 
mentioned previously, we assume a flat Universe with Q° M = 0.3 and Q° D = 0.7. For the 
sake of illustration, we split the DE density into its individual components as follows: 
= 0.65 and fl x = 0.05 (which, of course, add up to the value of fi° D ). In Fig.|3|, we 
show the regions of (wx, v) plane that fulfill the partial constraints: F < 0.1 (Fig. |^a) 
and F > —0.1 (Fig. ^b), as well as the more restrictive region \F\ < 0.1 (Fig. |3]c), 
which represents their intersection. All points in these figures automatically satisfy the 
nucleosynthesis bound ( f4.12|) . This means that, in this case, the F-test is already more 
restrictive than the nucleosynthesis bound. 

We note that wx is less constrained when v = 0, i.e., when the AXCDM model 
behaves as a quintessence (or phantom) model with a cosmological constant - see 
( f4.7| ). In this case, we have —14 < wx < —0.3. Similarly, when wx is very close to 
the value wx = — 1, the cosmon behaves as a cosmological constant and the range of v 
is almost unconstrained by the F-test. This scenario effectively corresponds to having 
a strict cosmological constant (in this case, the cosmon), together with a variable 
cosmological term, A. These unconstrained situations will change dramatically when 
the other restrictions (particularly that of non-cosmic coincidence) are also imposed 
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Figure 3: F-test condition ( [2.23D for the AXCDM model assuming flat geometry with = 
0.3 and individual DE densities = 0.65 and fl x = 0.05. Shown are: a) region of the (z/, wx) 
plane in which F < 0.1; b) the corresponding region where F > —0.1; and c) the allowed 
region by the F-test (|F| < 0.1), i.e. the intersection of regions a) and b) . 



in combination with the F-test. We have already seen, in the previous section, that 
wx = — 1 is actually forbidden. However, the F-test, alone, (as a strategy to determine 
the principal restrictions due to structure formation) is able to highly constrain regions, 
which the other conditions are not able to do, as we will see below. 

By assuming that Q° M = 0.3 for the present matter content of the Universe (which 
can be deduced from LSS observations alone) and that the Universe is flat, the AXCDM 
model is left with three free parameters, namely, wx, v and Q A . The cosmon density 
Q x is, then, no longer independent, since it is fixed by the cosmic sum rule ( |4.6|) . In 



a previous work MM , the parameter space of the AXCDM model was constrained by 
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imposing the three conditions ( |4.9| ), ( |4.1U| ), and ( |4.12| ), emerging from three important 
physical requirements: 



Condition 1: Nucleosynthesis bound. 

The expansion rate depends directly on the amount of DE. If we are not to spoil 
the standard Big Bang predictions (about, e.g., light element abundances), the 
DE density at the nucleosynthesis time should not be very large. Specifically, it 
is required that p7| 



&dM = < 0.1 , (5.4) 

1 + r N 



which is essentially equivalent to 

' Pd{z) 



<0.1, (5.5) 

(z=z N ) 



Pm(z) 

where zn ~ 10 9 is the cosmological redshift at the nucleosynthesis era. 
Condition 2: Stopping condition. 

As commented in section f|, in the AXCDM model, the ratio r(a) of the DE 
to matter density may exhibit a maximum. This feature is related to a future 
stopping and subsequent reversal of the expansion of the Universe, expressed by 

(EH)- 

Condition 3: Low maximum of the ratio r(a). 

For a solution, or at least a substantial alleviation, of the coincidence problem, we 
must further require that the stopping point of the expansion of the Universe is 
preceded by a sufficiently small maximum value of the ratio r(a), defined in ( |2.15| ). 
In the standard ACDM model, r(a) is unbounded and its present value, r ~ 1, is 
related to the recent transition from decelerated to accelerated expansion. Thus, 
within the standard ACDM, we can conclude that we are living in a very special 
epoch, i.e., very close to the transition epoch. Alternatively, if r(a) remains 
bounded and sufficiently small for the entire history of the Universe, as can be 
the case in the AXCDM model, the fact that r$ ~ 1 should no longer be regarded 
as a coincidence since the relation r ~ 1 could hold for most of the lifetime of 
the Universe. Therefore, in order to accommodate this appealing feature in the 
AXCDM model, we search for points in the parameter space that not only allow 
for the existence of a maximum for r(a), but also those points, for which the value 
of the maximum is sufficiently small. Specifically, we express quantitatively this 
condition by (f4.9|). 
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Conditions 2 and 3 are actually related since the existence of a maximum for r(a) 
is correlated with the existence of stopping. We have, nevertheless, kept them 
separated in order to stress that the stopping condition is not sufficient to solve 
the cosmic coincidence problem, it is only a necessary condition. We still have to 
demand a "qualified stopping", i.e., one preceded by a sufficiently low maximum. 
Conditions 2 and 3 could then be unified and be collectively referred to as The 
ability to solve the cosmic coincidence problem. 

The subset of points satisfying these three conditions was found in [Jy| to constitute a 
volume in the 3D parameter space, shown in Fig. ^ (upper panel). In this paper, we 
have incorporated a fourth and very important condition in the numerical analysis: 

• Condition 4-' Consistency with the data on structure formation. 

We have already explained in detail, throughout this work, the way in which we 
have included this condition in our analysis, namely, through the effective EOS 
approach and the implementation of the F-test, ( |2.23| ). This has been one of the 
principal aims of this work. 

There is yet one more observational requirement that can be demanded: 

• Condition 5: EOS condition at z = 

The value of the effective EOS parameter should behave as a CC in the recent 
past. However, the effective EOS of the AXCDM model (|5.1|) does not automat- 



ically satisfy the "CC condition", w e (z) — > —1 for z — > 0, as does the effective 
EOS of the running ACDM model (|3.7| ). Therefore, we wish to make sure that 
this condition is indeed satisfied by the AXCDM model, within the limits of the 
latest observational data. We normalize this EOS bound at z = 0, which is 
elaborated below. The quantitative restriction associated with this condition is 

\l + u e (z = 0)| < 0.3. (5.6) 



This can be justified from recent experiments, which strongly suggest that the EOS 
parameter should be close to —1. For instance, the combination of WMAP and the 
Supernova Legacy Survey (SNLS) |36j data (under the assumption of spatial flatness) 
yields §: 

u e = -0.9671°;™- (5.7) 
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Figure 4: Upper panel: Volume of points allowed by Conditions 1-3 described in the text, 
i.e. nucleosynthesis plus the two restrictions associated with the solution of the coincidence 
problem; Lower panel: Volume of points allowed by all constraints discussed in the text, 
namely Conditions 1-5, which include the LSS constraint and the EOS restriction at z = 0. 
We assume flat space geometry with Q° M = 0.3 and Q^, = 0.7. 

Even without the prior that the Universe is flat, the value of the EOS parameter 
preferred by WMAP, large-scale structure and supernovae data is still very close to 
that of a cosmological constant || , 



These observational results assume that the EOS parameter is constant and, therefore, 
they are not directly applicable to the AXCDM model. Nevertheless, we prefer to 
adopt a conservative point of view and impose the additional constraint from (|5.6|) , 
which we referred to as condition 5, to our model. 

Remarkably enough, after imposing the five conditions listed above, we are still 



uj e = -1.08 ±0.12. 
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Figure 5: z^-slices of the physical 3D region of the AXCDM parameter space (Fig.Q, lower 
panel) for v = i/q, 3z^o and — 2vq. Points in the solid-shaded area fulfil the nucleosynthesis 
bound and the two conditions associated with the solution of the coincidence problem ( Con- 
ditions 1,2 and 3 in the text). The dotted region consists of points allowed by the LSS data 
{Condition 4)- Points inside the dashed lines satisfy Condition 5: |l+o; e (0)| < 0.3 (left panel) 
and |1 + w e (0)| < 0.15 (right panel). The physical region (the darker one) in each case is the 
common overlap. Note that quintessence-like cosmons are forbidden if |1 +u; e (0)| < 0.15. 
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left with a non-empty volume of points in 3D parameter space that satisfy all of them. 
This can be seen in Fig. |] (lower panel). Comparing with the upper panel of this figure, 
we see that the new constraints greatly diminish the final volume of points allowed in 
the parameter space. This means that conditions 4 and 5 are very restrictive when 
combined with the first three conditions. 

If we wish to better visualize the effect of the different constraints, it is more 
convenient to study two-dimensional slices of the final 3D volume in Fig. ^ (lower 
panel). In order to do this, we fix one of the three free parameters. In Fig. |5], we plot 
three different //-slices of the final volume: v = z/ , 3z/ , — 2u Q , where z/ is the canonical 
value, defined in (|3.3j). In Fig. || we plot three 0° -slices: Q A = —0.65, 0.8, 1.2. Let 
us discuss the //-slices and restrict ourselves, for the moment, to the panel on the left 
of Fig. |5]. We see that both signs of v are possible and that the cosmon may be both 
quintessence-like [wx ^ — 1) and phantom-like [wx < — 1), although the former is only 
possible for positive v, according to Fig. || (middle figure on the left). This feature does 
not depend on this particular i/-slice [y = 3z/o), it is general. It can be confirmed from 
Fig. |7|, where the projections of the bulk volume of final points (lower panel of Fig. |j) 
onto the three possible planes are shown. Also general is the property that phantom- 
like cosmons are compatible with positive (greater than +0.7) values of Q° A , whereas 
quintessence-like cosmons demand a negative energy density of the vacuum (cf. Fig. |7|). 

Since the precision of observational data will increase in the future, we have also 
studied what happens if we make the condition on u; e (0) more stringent than condition 
(|5.6|), namely, 

|1 +w e (0)| < 0.15. (5.9) 

The result of applying this tighter constraint is that the volume of points allowed 
becomes further reduced, as can be seen in the panel on the right of Fig. [|. The most 
conspicuous effect is that, now, the X component can only be phantom-like [wx % — 1); 
quintessence-like cosmons (wx > — 1) are no longer permitted. 

Consider now the f2°-slices (Q A = —0.65, 0.8, 1.2) of the final volume in Fig. ^. In 
general, only a small subset of all the points allowed by the F-test is allowed by the 
full set of constraints. Conversely, the area allowed by the old constraints ( Conditions 
1, 2 and 3) is highly restricted by the F-test. For example, the Q° A = —0.65 slice shows 
a critically constrained case, where the concurrence of the five conditions, particularly 
the EOS restriction at z = (Condition 5), leave a very small area allowed in that 
slice. 

Let us now consider the projections of the final allowed volume in Fig. |7|. The signs 
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Figure 6: Three fi^-slices of the physical 3D region of the AXCDM parameter space for 
ftO = -0.65, 0.8, 1.2. The shaded region is allowed by Conditions 1-3 in the text. The dotted 
one is allowed by Condition 4 (the F-test), and the region between the vertical dashed lines 
is allowed by Condition 5. The = 0.8 plot is not restricted by this condition. The darker 
areas in the three plots are the final allowed regions. 




of v and are correlated, as shown in Fig. [7|b. Negative values of are compatible 
only with positive values of v. Similarly, positive values of admit mostly negative 
values of z/, although there is still a range of v > allowed values, which are of order 
of the canonical value uo ~ 10 -2 . From Fig. [7|c, we see that the signs of 1 + wx and 
fl\ are also correlated. Thus, quintessence-like cosmons (1 + wx > 0) demand that 
< 0, whereas phantom-like ones (1 + wx < 0) require that Q° A > +0.7. However, 



it should be stressed that if we would apply the tighter condition, indicated in (|5.9|) , 
only the phantom-like character is possible for the cosmon; in other words, the thin 
quintessence regions in Fig. 0a and Fig. 0c then disappear altogether. 

The reduction of the original region, permitted by the first three conditions, to the 
final region, filtered by all the five conditions, is best seen, comparing the projection 
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Figure 7: Projections of the final 3D physical region of the AXCDM model (Fig.|j, lower 
panel) onto the three different planes. The vertical and horizontal dashed lines mark the 
values v = 0, wx = — 1 and Q A = +0.7 wherever any of these parameters appear. 



plots in Fig. |7| with the old projection plots in Fig. 3 of ]18 |. For instance, the present 
range of values for the v parameter become bounded from above and below and it is 
no longer possible to have \v\ > 0.5. Another dramatic restriction that occurs in the 
present analysis is the near exclusion of the region where the cosmon is quintessence- 
like, the very small region of the allowed parameter space where wx ^ — 1- In this 
region, v cannot be negative or zero (or even close to zero) and must be v > 0.1. 
This is clearly seen in the projection plot of Fig. |7|a. We conclude that in most of the 
parameter space, the cosmon X behaves phantom-like (wx < —1)- Since, however, 
Q° A > +0.7 in this region and the current value of the total DE density is fixed at 
= +0.7, the cosmic sum rule Q4.6Q indicates that in most of the allowed parameter 
space, Q° x < 0. Put another way, the cosmon behaves mostly as phantom matter. 
This was to be expected since phantom matter satisfies, as noted previously, the strong 
energy condition ||18|| , and therefore it effectively behaves like additional matter, helping 
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Figure 8: The effective EOS ( O ) of the AXCDM model as a function of the redshift for a 
set of allowed values, see Fig.0: (a) wx = —1-6, = 0-80 (phantom-like cosmon) for three 
values of v\ (b) wx = —0.8, Q\ = —0.3 (quintessence-like cosmon) for three positive values 
of v (no negative values allowed in this case) 

to stop the expansion of the Universe at some time in the future. 

We emphasize that this preference for the phantom matter character of the cos- 
mon entity does not necessarily translate to the overall effective EOS of the AXCDM 
model. This has already been discussed after Eq. (|5.3| ). The quintessence- like beha- 
vior [dVlo/da < 0) or phantom-like behavior (dQo/da > 0) of the total DE of the 
AXCDM model depends only on the sign of the product (1 + wx) fix m the regions 
where Qr> is positive. Let us consider some examples, particularly those in Fig. |]a, 
which correspond to allowed values of the parameters. Since, in this case, Q x = —0.1 
is small, from Eq. ( [5.2|) we can see that for sufficiently large z (small a), the last term 
on the r.h.s. is the dominant one (using the fact that e 1). Therefore, we have the 
following rule of thumb: the total DE of the AXCDM model behaves as quintessence 
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(w e > —1) or phantom {w e < —1), if the sign of v (1 +wx) /wx is positive or negative, 
respectively: 

l + ^^O^ ^ + ^^ Q. (5.10) 

w x 

For example, if X behaves as phantom matter (which is the preferred state of the 
cosmon, according to our analysis), the AXCDM model can be phantom-like if v < 0, 
but behaves like quintessence if v > (cf. Fig. §a). When the first term on the r.h.s. 



of Eq. (ET2) becomes important, the rule applies only if z is sufficiently large, as it is 
the case of the examples in Fig. ^b, where fl x = 1. These examples illustrate that 
the overall behavior of the effective EOS of the AXCDM model is not dictated by the 
cosmon EOS alone, but depends on other parameters, in particular v. 

Finally, we may compare Fig. |8] of the AXCDM model with Fig. [I] of the running 
ACDM model, studied in the previous section. For the AXCDM, the allowed values 
of v are larger. As a result, the departure of the effective EOS of this model with 
respect to the CC divide, w e = —1, can be significantly greater in the redshift ranges 
that are accessible to the next generation of supernovae experiments. It is, therefore, 
possible to investigate deviations of w e from -1, predicted by the AXCDM model, in 



these experiments, particularly by SNAP p5 |. For instance, in Fig. ga we see that for 
v = —pq, we have the predicted value w e (z = 1.7) = —1.06 and for v = +vq, we have 
w e (z = 1.7) = —0.92, corresponding to the highest redshift accessible to SNAP. The 
departure from —1 is significant large to be detected. At z = 2, which will be accessible 
in the future, we have w e = —1.10 and w e = —0.90 for the same values of v. 



6. Conclusions 

In this paper, we have explored the growth of matter density perturbations for two run- 
ning A models in the literature, the running ACDM model [|7|, |8], 11] and the AXCDM 



model [TE|, pU[ . These models offer alternative explanations for the dynamical dark 
energy (DE), beyond the usual proposals based on quintessence ideas. In particular, 
the AXCDM model constitutes a promising cosmological framework, of a very general 
nature, that also has the capacity to try to understand the conspicuous cosmic coin- 
cidence problem, namely the perplexing coincidence of finding ourselves, at present, 
in an expanding Universe, where the amount of dark energy turns out to be of the 
same order as that of matter. The cosmic coincidence is a riddle, wrapped in the 
polyhedric mystery of the Cosmological Constant Problem [|I|] , which has many faces. 
In this case, the conundrum is to understand the following situation. The density of 
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matter is continuously decreasing with the expansion of the Universe, whereas the DE 
energy remains constant in the standard ACDM model. It is totally inexplicable to 
understand from this model, the fact that the cosmological constant line crosses the 
ever-falling matter density curve, precisely now. The AXCDM model may be able to 
provide a clue to the resolution of this enigma, due to the dynamical interplay of the 
cosmon X and the variable A, which together, form a composite dark energy medium. 
This dynamics insures that the ratio of the DE to the matter density stays bounded 
and that its value, at present, is not very different from the value it will have in, say, 
the next Hubble time. 

An important aspect of the AXCDM model that has been yet to be investigated 
is its consistency with the LSS data. In this paper, we have undertaken a thorough 
study, in an attempt to try to answer this crucial question, namely, is there a non- 
empty region in the 3D parameter space of the AXCDM model, capable of solving the 
cosmic coincidence problem and, at the same time, being consistent with the known 
data on structure formation? Our study shows that the answer to this question was, 
in fact, positive. 

The "effective approach" that we have used here is based on three essential ingre- 
dients: i) the use of the effective equation of state (EOS) representation of cosmologies 



with variable cosmological parameters [|T^, [13| ; ii) the calculation of the growth of mat- 



ter density fluctuations, using the EOS of the DE |28|, p9| ; and iii) the application of 
the "F-test" |30|, |31[ to compare the model with the LSS data. This three-step method- 
ology turned out to be a streamlined strategy. Even if it is not a perfect procedure 
to estimate the restrictions that structure formation imposes on a given model of dark 
energy, it is nevertheless an economical and efficient method to encapsulate the essen- 
tial findings of the full-fledged approach. That this is so can be argued on the grounds 
of the various existing cross-checks on the constraints that LSS data impose e.g. on 
the running ACDM model, the first dynamical A model that we addressed in this 
study. Our "effective approach" provides a noticeable consistency with the complete 
calculation of matter and DE perturbations presented for the running ACDM model 
in Ref. [2tJ. We are, therefore, confident that the same procedure is able to capture 
the main restrictions that the present data put on the parameter space of the AXCDM 
model. 

In view of the consistency between the solution to the cosmic coincidence problem 
proposed by the AXCDM model and the data on structure formation, this model is 
substantially reinforced. It emerges as a very strong candidate for a possible solution 
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to the cosmic coincidence problem. The physical region of its parameter space turns 
out to be compatible with all cosmo logical data known at present. Furthermore, we 
have shown (cf. Fig. |J) that the model predicts non-trivial observable features in the 
EOS of the dark energy. Most important, these features can be accessible to the next 



generation of supernovae experiments, like DES and SNAP |35|]. This model, therefore, 
has the ability to solve some of the important problems of modern cosmology and, in 
addition, makes predictions which can be tested by observations from experiments just 
around the corner. We eagerly await the possibility of confronting the AXCDM model 
with these observations. 
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